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^ , Abstract 

^ . Using a set of new odderon states, we calculate their contribution to the difFractive 77^ photo- 

' and electroproduction process. Compared to previous simple 3-gluon exchange calculations we 

. find an enhancement of about one order of magnitude in the cross section. The i-dependence of 

fSJ ' the cross section exhibits a dip structure in the small t region. 

o : 
^ ■ 

O ; 1 Introduction 

^ : 

Q-i! The existence of the Odderon Q , the partner of the Pomeron which is odd under charge conjugation 

Qh' C, is an important prediction of perturbative QCD. In the leading order, the Odderon appears as 

^ ■ a bound state of three reggeized gluons. Its experimental observation is a strong challenge for the 

• • , experimentalists. Promising scattering processes where the exchange of the Odderon may be seen 

. ^ I include the difference between the p — p and p — p cross-sections and the diffractive production of 

^ ■ particles with a C-odd exchange, such as photo- and electroproduction of pseudoscalar mesons (PS). 

From the theoretical point of view, some of the latter processes are of particular interest, namely 
those where the presence of a large momentum scale provides some justifications for the use of 
perturbative QCD. This includes, in particular, the diffractive production of charmed pseudoscalar 
mesons, for example the ijc- Correspondingly, a large amount of literature has been devoted to this 
class of diffractive processes. For large photon virtualities Q^, for heavy mass PS mesons (such as 
rjc), and for large momentum transfers the relevant impact factors for the transition 7(7*) — >PS 
have been be calculated perturbatively Q. As to the Odderon exchange, first studies have used 
the simplest form, the exchange of three noninteracting gluons in a C=-l state. Another line of 
Odderon investigations, persued by the Heidelberg group, uses a non-perturbative model for the 
Odderon, based on the idea of a "stochastic QCD vacuum" Numerical estimates for the cross 
sections turn out to be very different in these two approaches: the nonperturbative Odderon models 
tend to give substantially larger cross sections Evidently, to guide the experimental search 

for signatures of the Odderon we have to clarify these discrepancies. In this paper we follow the 
perturbative approach and make use of the recently discoverd new Odderon solution. 
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The perturbative QCD Odderon has a rather long history. After several variational studies ||5|,|^ 
a first analytic solution to the Odderon equation was constructed by Janik and Wosiek and 
verified by Braun et al |8|. It belongs to the lowest non-zero eigenvalue of the conformal integral 
of motion, Q^. This solution has an intercept slightly below unity and, most important, vanishes if 
two of the three gluons are at the same point. This property leads to the disappointing result that 
this solution cannot couple to the perturbative 7 — >PS vertex and therefore is irrelevant for photo 
and electroproduction of PS mesons. It may, however, play its role in purely hadronic processes, 
such as pp or pp scattering. Recently a new solution for a bound state of three reggeized gluons 
with the Odderon quantum numbers has been found |^], which is quite different from the previous 
one. It corresponds to Qa = 0, has intercept unity, and, most important, it does couple to the 
7 — >PS transition vertex. The structure of the wave function of this solution is rather peculiar, so 
that one may expect substantial changes compared to the exchange of three noninteracting gluons. 
The purpose of this paper is to perform an analytic and numerical study of the exchange of this new 
odderon solution for the diffractive ?7c photo- and electroproduction. We follow the approach of 
by replacing the three noninteracting gluons by the new Odderon state. As the main results, our 
cross sections are an order of magnitude larger than those of QQ; also, the t-dependence of our 
cross section exhibits an interesting dip structure which is not present in the case of noninteracting 
gluons. 



2 The perturbative QCD Odderon which couples to 7-PS transi- 
tion. 

Let us first briefly recapitulate the main properties of the new odderon solution. For details we refer 
to ||9|. We remind that this solution has been shown |]l^] to be a particular case of a more general 
class of solutions: in the large A'c-limit there exist relations between eigenstates of different number 
of reggeized gluons, which connect solutions of different symmetry properties, and the Odderon 
solution is the simplest case that satisfies these relations. 

In momentum space the Odderon wave function is constructed from the known Pomeron solu- 
tions [12 1, which have the eigenvalues 

X{v,n) = as ^2^(1) - i^ C'^}'^^ - "L^"'^ " i^)] 1 = (1) 



One proves that 



(123) '^l'^2 



(2) 



indeed satisfies the Odderon equation and has the same intercept (|T]), provided the Pomeron wave 
function is odd under the interchange of its arguments. This restricts the values of n to odd 
numbers. The normalization factor c in (|2|) can be chosen in such a way that the Odderon wave 
function will have the same norm as the Pomeron function 

= ^^{u;n)^^{u',n')^ ^ w{u,n)5{u - u')5nn' , (3) 

where the w{v,n) are known |11|. In (^) the scalar product is defined as the integral of the two 
wave functions in momentum space, where the bra- vector has to be amputated (marked by a tilda). 



2 



Condition (^) leads to 



When including in (^) the colour structure dabc one should change the normalization by a factor 
VK/V{N--A){N--1). 

In order to construct the full Greeen function we clearly need to know the complete set of 
solutions of the Odderon equation. At the moment we only have the symmetric solutions of Q 
and the new solutions (^). The former ones are orthogonal to the photon impact factor and so 
irrelevant for our problem. At present we do not know if any other solution exists, apart from 
(|2|), which couples to the 7 — PS transition vertex. So our results are strictly speaking restricted 
to the contribution of the exchange of the Odderon states with the wave function (refoddwave) . 
Normalizing the Green function to reduce in the small coupling limit {as 0) to 

and having in mind @ we find the part of the Green function corresponding to (^) in the form: 

G3{y\ki,k2,k3\k[,k2,k'r^) = 

omJ-°- [u^ + (n-lY/A\[v^ + {n + lf/Ay ^^^i' ^^2, ^3)^ [k„k2,k,). (5j 

3 The BFKL function in the momentum space. 

In this section we present the BFKL Pomeron eigenstates in the momentum representation, which 
we need in order to construct the Green's function in (|5|). This wave function is well known in the 
coordinate space |[T^ where its form is dictated by conformal invariance 

Ei''^^\r,o,r,o)=(^y(^)\ (6) 
\rior2o/ \rior2o/ 

where rio = ri — etc, h = (1 + n)/2 + iz^, h = (1 — n)/2 + ii^, and the standard complex notation 
for the two-dimensional vector is used on the right-hand side. Fourier transforming to momentum 
space one finds (see Appendix 1) 

E^nik.M) = E^j^ik,, k2) + El^[k,M), (7) 

where the first term denotes the analytic contribution, and the second one stands for the 5-function 
terms. The analytic part is given by 

Ef^j^ikiM) = c(x{kiM) + {-iTx{k2M)). (8) 



where h = {\ + n) /2 + iv and h = {\ — n) /2 + iv are the conformal weights. The coefficient C is 
given by 



C =^^j^hh{l-h){l-h)T{l-h)T{l-h). (9) 
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The expression for X in complex notation is given in terms of the hyper geometric functions 



h~2 



The (5-function part is 



F {l-h,2-h;2; 



F [l-h,2-h;2; 



(10) 



^ 9 = fci + fe2. (11) 



Let us note that, when the Pomeron couples to an impact factor of a colourless object, (5-function 
terms in the pomeron wave function do not give any contribution. However in our calculations it 
turns out that these terms play an important role. 

In order to calculate the couplings of the exchanged Odderon to the proton and to the 7 — > ?/c 
impact factors we need to know the Pomeron wave function E'^'^'") in momentum space for n = ±1 
and around the value u = 0. We start with the analytic part; we need to look at (10) which, for 
n = 1, reads 



X{ki,k2) 



iu—l 



F —iv, 1 — iv] 2; 



F\l-iu,2- iu\2\ 



(12) 



and we perform a Taylor expansion around the point 1^ = 0. The complicated expression (|^) is 
drastically simplified at small values of u. In lowest order it is linear in z/: 



Et{ki,k, 



1 



2-11^ q \k1k2 ^2^1 



(13) 



This function is odd in the azimuthal angle (i.e. antisymmetric under 93 ip+ir). So it is orthogonal 
to the two impact factors which are azimuthally even. For this reason a non-zero contribution only 
comes from the terms quadratic in i'. Omitting those of them which have the same structure as 
dg) we find 



E^{kiM 



2 ^ 



27r2g 



1 



1 



1 



In fcf - In fc^ + To- - 7^ In g + 



k 



k 



k2kl kik\ 



qlxiq 



The (5-function part is finite for v ^ 



i 1 
2^'q' 



The constant C in (^) (for n = 1) becomes: 

c ' 



(47r)2 

To construct the Green's function one has also to consider 

Xii^,±l) = -iy^2aa3)+0{u^) 
and the i^-dependence in the denominators of the integration measure in (^): 

[1/2 + (n - 1)2/4] [1/2 + + 1)2/4] ^ ^2(1 ^ ^2) 



(14) 
(15) 

(16) 

(17) 
(18) 
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4 The transition amplitude 



We study the process of the diffractive photo- or electroproduction of rjc on the proton: 7*p rjcP- 
It will be assumed that the proton remains intact, although it would be rather easy to include also 
its low-lying excitations. The differential cross-section is given by the formula 



,,(7(7-)+P-,. + P) = ^^i:i-4f, (19) 

1=1 

where A^, i = 1,2 is the photoproduction amplitude for a given transverse polarization i of the 



photon i = 1, 2. The electroproduction cross-section can be obtained from (19) in a trivial manner 
(see @). The photoproduction amplitude is given by a convolution of the two impact factors, 
$p and for the proton and for the 7 — > 7?c transition, resp., with the Odderon Green function: 

"511 ($;|G3|$p). (20) 



32 6 3!(27r)8' ^ 

We shall assume that the cm. energy squared, s, is much greater than the scales in the transition 
vertex: s » Q'^,m'^,t. Using the definition of [Q], namely a; = {m'^^+Q'^)/{s+Q'^) and y = logl/x, 
we are in the low-x limit x ^ 1. 

The matrix element on the rhs of (^0|) involves the integration over v coming from the Green 
function (|5|). This integration can be done in the saddle point approximation, since y >> 1. The 
leading contribution will obviously come from the smallest values of |n|, i.e. n = ±1 and small u, 
when the integral acquires a form 

J due-'^'fyiiu), /5 = 2a,C(3) (21) 

and I(i^) denotes the nonexponential part of the integrand. Since we expect the dominant contribu- 
tion to the cross section to come from the kinematical region where t is not large, there is one more 
large momentum scale in the problem, apart from the energy ^/s, namely, M = \/Q''^ + 4m^, which 
provides some basis for the use of perturbation theory. In principle, the position of the saddle point 
i^s depends upon the relation between these two large scales. However in our kinematical region of 
small X, M is much smaller than ^/s, so that I'g is close to zero. As a result, we have to calculate 
the matrix elements of the two impact factors with the Odderon wave function with |n| = 1 and 
small u. 

The impact factor corresponding to the transition ^ ^ rjc can be calculated perturbatively. 
Referring the reader to the original paper |Q] for the details, we only quote the final result: 

^^=b,..^(\- (fci + fc2 - fc3) • g ^ \ .22) 

^ \(^)Q' + 4m2 + (/ci + fc2-fc3)2 Q^ + Aml + q^j- ^ ' 

with q = ki + k2 + k^ and 

16 ol , 

0= —ecgs^rn^JjQ. (23) 
vr 2 

Here Cc = (2/3)e is the electric charge of the charmed quark, and Qs is the strong coupling constant. 
The constant 60 can be determined from the known radiative width r(r/c — > 77)= 7 KeV: 



IGtt'^ vrT 
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Following we denote 

, _ ik-k')ik + k') 

Using the form (||) of the Odderon wave function one can show that the convolution of the Odderon 
wave function with the impact factor for the transition 7 — > ?/c reduces to a convolution of the 
function ip and the Pomeron function E 

(f;^,J*(^'")) =5ei,4-l- / cPk^{k,q-k)E^''^^\k,q-k) 
' ' civ, n) J 

where the Pomeron function is given in (0), and we have rescaled the momenta in the integrand by 

Note that ( p6|) is sensitive to the (^-function term present in the Pomeron wave function. In 
fact in ( p^ ) has been proven by explicitly using the Pomeron wave function in the coordinate 
representation, which does contain such terms. 

We could have chosen to calculate in a different way, directly from (^): 

($;^,J^(^'")) =3c J dHE*(''^^\l,q-l)g{l), (27) 
where function g{l) is defined by 

gil) = J d'k l-k,q-l), (28) 

factor 3 comes from the symmetry of the impact factor and ki = k, k2 = I — k and fca = q — I. 

Now the (^-function terms in E give no contribution. Indeed, g{0) = because of the factor 
present in (p8[), and g{q) = as the impact factor vanishes for k^ = 0. As a result, in the calculation 



of (P^), one may freely add or subtract such 5-function pieces in the Pomeron wave function |14|. 
At V = ±1 and small i' factor c behavies as l/v (see (^) and (^) ). Therefore, to find the 7 — PS 
form-factor in the lowest (first) order in u one needs to know the BFKL function E^'^"^^^ in the 



second order if one uses (|27D , but only in the zeroth order if one uses (26). One can show that 
both ways lead to the same answer. However, ( ^7|) needs numerical computation, and (^) gives 
the result in a trivial manner due to the simple structure of the zeroth order BFKL function (|l5|) . 
At |n| = 1 and ^ we find 

y(o,±i) M ^ _ VWM^ (29) 
In Fig. 1 we show a plot of this function in the region < |t|/M2 < 10. We have verified by 



numerical computation that, starting from (27) and using for E*^'^'"'^ only the analytic part, E2, 
we get exactly the same answer for ^ 0. 

The proton impact factor is non-perturbative. We use the parametrization proposed in 

3 



F{q, 0, 0) - ^ F{ki, q - k^,0) + 2F{ki,k2, fcs) , (30) 



i=l 



6 



^(0,±1) 
7 

0.2 



0.1 
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Figure 1: Numerical results for the coupling of the Odderon to the 7* 
in (p6[), as a function of the scaled variable x = |t|/(Q^ + 4m^). 



r/c impact factor (defined 



with 



F{ki,k2,k3) 



2a' 



k2? + (fc2 - ^3)2 + (fc3 - k^] 



(31) 



2a? + (fci 

d = 8(27r)^5^, and the scale parameter a = mp/2. From the comparison with the two gluon 
exchange model for hadronic cross-sections the authors of |||] estimate ^■^/(47r) = 1. The impact 
factor ( |30[ ) satisfies the basic requirement that it vanishes when any of the three gluon momenta 
goes to zero. The calculation of the scalar product of the Odderon wave function with the proton 
impact factor is more cumbersome, since it amounts to the integration over the three-gluon phase 
space. Using the symmetry of the proton impact factor and the explicit form of the Odderon wave 
function (0) one finds 



=3cd j d'^lE*^'''''\l,q-l)f{l) =cd 



1 



(2a2)i/2 



. [if 



2\* 



y{u,n) ( ^ 



where 



d'k- 



k\i-kY 



[F(q,0,0) 



3 

^ F(kj,q - kj,0) + 2F{ki,k2, k^] 



(32) 



(33) 



and, as before, fci = k, k2 = l — k and k^ = q—l. To obtain the last expression in ( [3^ ) we have again 
rescaled the momenta in the integral with respect to q and extracted the factor ii/2/27r2 having 
in mind that at small z/ the contribution to ( |3^ ) starts with the second order term of the BFKL 
function (0). Note that the integral (33) is infrared finite, since the square bracket vanishes if any 
of the gluon momenta go to zero. However, individual terms inside the square bracket are infrared 
divergent. Clearly, as in (p7|), only the analytic part of the Pomeron function contributes to 



Let us consider the rhs of (^). Two of the five terms are simple, since the -F's do not depend 
on the integration variable. They give 



/i = (F(q, 0, 0) - F{q - I, Z, 0)) J d^k-^ 



I' 



F(q,0,0) - F{q -1,1, 0))2TTln 



The two remaining terms in the sum inside the square brackets in (33) give identical contributions. 
Their sum is given by 



/2 



I' 



1 



--fl2 / d'^k- 

3 J ''k\l-k)Hk-q/2y + e' 



(35) 
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Figure 2: Numerical results for the coupling of the Odderon to the proton, defined in (^2|), as a 
function of the scaled variable x = \t\/2a?'. 



where e = + The last term in the sum gives: 



k\l-k)Hk-l/2f + 5' 



where 5 = I — Iq + {q + a )/3. Both integrals, and (pq), thus reduce to a generic one 



d'k- 



1 



k^ + m?)[{k - If + m^][{k - pf + 62] ' 



(36) 



(37) 



corresponding to a general two-dimensional triangle diagram. We have introduce here a mass m 
as an infrared regulator. Using the standard Feynman parametrization this integral can be trans- 
formed into an one-dimensional one which, after separating the infrared divergent contributions, 
can be done numerically. Some details about this procedure are discussed in Appendix 2 (note that 
(|3^ ) actually corresponds to a collinear case which, in principle, can be calculated analytically). 

The resulting function /(Z) = fi + f2 + is used to calculate the integral (^). Here only the 
analytic part, £'2^, contributes. The calculation has been done numerically. We present the results 
for vjP'^'^'' in Fig. 2. One observes that Vp changes sign at \t\ « 0.07 GeV^. As a consequence, the 
cross section will vanish at this point. Of course, this property is literally true only in the limit of 
asymptotically large energies where one can neglect the contributions of all other states v ^ Q and 
Inl > 1. 



5 Numerical Results and Discussion 

To find the final cross-sections from ( [l9|) and (|20|) , we do the saddle point integration over v, then 
take the square module of the amplitude ( [20| ) and do the sum over the polarizations in (|l^). The 
latter step provides a factor l/|f| (from the prefactors in (|22|)). The normalization factors c(z^, n) 
of the Odderon solution which are contained in the Green's function (^) cancel when, in (p6|), the 
scalar product of the Odderon wave function and the photon impact factor is computed. Collecting 
all pieces of our cross section formula we find 
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The differential cross sections for the two cases = and 25 GeV^ and ^/s ~ 300 GeV are 
shown in Fig.3. We have taken at the scale + (in at = the scale was m^). 
Independently of the cross sections show a dip at small |t| ~ 0.07 GeV and a maximum at 
\t\ ~ 0.22 GeV (roughly of the order of (mp/2)^. The dip comes from the zero present in the 
Odderon-proton coupling (Fig. 2). Its origin seems to be related to the symmetry properties of the 
Odderon solution: as it can be seen from (|2|), the Odderon wave function is a sum of three terms, 
each of which contains an antisymmetric Pomeron eigenfunction. A similar structure is present in 
the photon impact factor (^), whereas the proton impact factor (^) is completely symmetric. The 
convolution of the Odderon wave function with the photon impact factor has no zero in t, whereas 
the convolution with the proton leads to such a zero. Whether this feature is an artifact of the 
simple model for the proton impact factor that we have used, or whether it represents a general 
property of the Odderon-proton coupling we do not know. In our calculation, this dip is present in 
the leading high energy approximation; it may be that at finite energies (e.g. at HERA) the dip is 
(partially) filled by the exchange of nonleading Odderon states. 

At t = the cross section vanishes, as in the case of a simple three gluon exchange. As one can 
see from Fig. 3, this happens at quite small t and cannot be seen in the figure. A more detailed 
analysis of the behaviour in the region of very small t requires, probably, a more accurate evaluation 
of the j/-integral in the 3-gluon Green's function. It would, however, not affect too much the dip 
structure or the value of the integrated cross section. 

1000 
100 

da 10 
1 



0.1 
0.01 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 I |^| 

Figure 3: The differential cross sections (in pb / GeV^). The upper curve refers to = 0. 

For the integrated cross sections we find 50 pb and 1.3 pb at = and 25 GeV^, respectively. 
Compared to the value of 11 pb predicted for = with a simple three- gluon exchange Q we find 
an enhancement of about 5 times. For = 25 GeV^ the differential cross sections in Q seem to 
indicate that one can simply scale the = cross section by 0.01 and obtain O.lpb. This implies 
that our cross-section is an order of magnitude larger. Note that compared to the simple three gluon 
exchange, we have a (weak) logarithmic suppression with energy. So the obtained enhancement 
effect is totally due to the coupling of our Odderon wave function to the impact factors. 

In conclusion, by comparing the exchange of three noninteracting gluons with the exchange of 
the new odderon solution we find that the interaction between the exchanged gluons leads to a 
significant change in the scattering cross section. However, despite this improvement in accuracy, 
our numerical estimate of the rjc cross section still suffers from a few theoretical uncertainties, in 
particular due to the Odderon-proton coupling. We believe that both the structure of the vertex 
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and its overall normalization should be checked more carefully. As a possible strategy, one might 
study the exchange of the Odderon in pp and pp scattering at large t where the use of perturbative 
QCD can be justified. Using the same model Q, a comparison with experimental data on the 
difference of pp and pp scattering fixes the overall normalization of the Odderon proton coupling. 
As to the general momentum structure of the vertex, the most sensitive test is the t-dependence in 
the small-t region: the presence of a dip would support both the structure of the Odderon-proton 
vertex and of the Odderon state used in our calculation. 
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Appendix 1. A momentum space representation of the Pomeron 
wave function in the non-forward direction 

The Fourier transform of the function E^'^'^^ given by Eq. (P) is defined by 



(27r)2 (27r)2 \r1r2 



Eu-HikiM) = ^77^4 eH---i+---^i . (39) 



1'2 



The integrand in ( p9| ) becomes constant as |t'i^2| 00 and therefore contains terms proportional 
to 5^(fci) or 52(^2)- Therefore one expects to find 

EhhikiM) = Ef:j^{kiM) + Eihiki, k2), (40) 

where the first term denotes the analytic contribution and the second the 5-like one. 

We shall at first compute the analytic part of the Fourier transform. Since we are dealing with 
a distribution, it is convenient to consider a new, regularized, object 

V/.aHa^'^l''^^) - J (2,^2 (2^)2 ^^^^^^H ^^^^^^^f ' (41) 

with an independent conformal weigth /13 for the ri2-terms. The integral (p|) is well defined for 
3f?(/i + h) <2, 5R(/i3 + /13) > —2 and 3ft(/i + h — — h^) > 0. Strictly speaking, the last inequality 
holds, for example, for the ^2 / case. When ^2 = one needs K(/i + h — hs — h^) > 2. In 
the calculation of the scalar product of the Pomeron function with test functions which vanishe at 
the points ki = or ^2 = we shall, therefore, consider the following prescription: {^\Ef^j^) = 
\im.fi^^fi{^\Efif^ fe-jfeg) = li™/i3^/i hghg)- This means that we shall be able to extract the 
analytic part of i^^- 

Introducing, for both external gluons of the Pomeron wave function, the complex variables 

hx iky 



2 



rx + iry = X + k 



iry = x-iy K = — 2-, (42) 

and remembering that h — h = n and h-^ — h^ = n^, we get a double integral in the complex plane: 

^Sh3(^l'^2) = / |^2|-h|^2|-h|^2^|h3,p^*n^*-n3^.(.in+.frJ+K..,+K^2) _ (43) 
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By using the integral representation 

1 



da a"-^e-"^ (3f?(x) > , > 0) (44) 



r{u) Jo 

we arrive at the representation 

/■OO 

Jo 

where, in order to fulfill the restrictions in Eq. (|4^), we need to require 

3ft(/i) > , 5R(/i3) < n > -1 , ng < -1 . (46) 

The next step is to perform a suitable change of variables in order to do the spatial integrations. 
Thanks to the holomorphic separability of the integrand, the natural choice would be to use rj and r| 
as independent integration variables. This can practically be achieved in the following way: first of 
all note that the whole integrand, regarded as a function of (xi, yi, X2, 2/2), is analytic. Therefore, we 
can analytically continue to complex values of, say, yi and 1/2, and rotate the respective integration 
paths by negative angles. This can be done by putting y = e~^^w : G 7?. for both j = 1, 2, and 
by simultaneously rotating the a- integrations in such a way that a = e*"^/? : P £ TZ. The generic 
exponential in the second line in (^5|) becomes 

which does not grow for large x and w provided < '7t/2 , \<j) — 29\ < tt/2. Choosing the extreme 
case (j) = 9 = tt/2 we set y = —iw , a = i(3 and perform the change of variables 

Pj = + ' Pj = - U = 1; 2) . (48) 

A further remark concerns the values of the complex momenta kj: the exponential in the first line 
in (E^) contains two factors 



J{ka:X+kyy) _ i{kxX—ikyw) _ ik^x kyW 



(49) 



In order to have a meaningful integral, the values of the y-components of the momenta must be 
imaginary. Hence we have to introduce a slightly different notation with respect to Eq. ( ^2| ) 

K= ^""'^^ ; r=:^2L±1^/k* ; K,RGn (50) 



and to consider k and k as independent variables. After this replacement we get 

47^(^1, K2;ki,K2) = ^r2YMrf-T)' rjP'"^^' (p-iP2)>-i2"^e^(^^^^+-^^-^)x 



/OO 
dpidp2 e-^(^l''l'^l+^2P2p-2+/33Pl2Pl2-^lPl-^2P2) _ (5^) 
-00 
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The last integral in the above expression provides two delta functions constraining 

Kl \ _ f Pi + (33 -/?3 \f Pl\ ^ f Pl\ _ 1 f (32+P3 P3 \f 1^1 



If, for simplicity, we restrict ourselves to positive values of the k's, it is apparent from the second 
equality in the above expression that the positivity of the /?'s forces the p-integrals to contribute 
only in the positive real half-plane. At this point we split the /^-integral into two pieces: the 
first takes into account the region p2 < pi and the second the remaining one p2 > pi- The first 
contribution reads 

dp J dp2{piP2TPu'e'^~^'''^~^''''^>< (52) 

d/3id/32d/33 {f3i(32f'^f33'''-H{[3ipi + ^^2 - i^iW2P2 - PsPu - '^2). 

We can factorize the above expression in the product of two independent integrals by means of the 
change of variables Aj = f3iPi : i = 1, 2, 3 , {ps = pu), which yields 

dpi / dp2 (piP2)""'^Pi2"'+'''e^^^''''+"''''^x 

^0 

/ dAidA2dA3 (AiA2)'^-U3 '^^^-^^(Ai + A3 - ki)6{X2 - A3 - K2) . (53) 
Jo 

The p-integrals are easily evaluated by setting x = P2/P1, which casts the inner integral into an 
integral representation of the confluent hypergeometric function iFi (see, e.g., Eq. (13.2.1) of [p!5| ). 
The outer integral then becomes a Laplace transform of the i-Fi times a power (see Eq. (7.621.4) 
of |l^). The result for the first factor of Eq. ( |53|) is (remember n — h = —h) 



^,3_,, r(l - h)ni + /.3)r(2 - 2h + /.3) 2H-H3-2 ^^^(1 -h,2-^ + h,-2-h + h3- -^) . (54) 

r(2-/i + /i3) Ki 



All the restrictions in the above formulas are automatically fulfilled because of Eq. (pq). The 
A-integral in Eq. (^) transforms into 

K^-'^-l^-l I' y-h-S-\l _ yf-\l + lLyf-\ (55) 
JO K'2 

which is just an integral representation for the hypergeometric function (see, e.g., Eq. (15.3.1) 
in [|l5|) and yields 

^ffl^.^-'^-l^-l ^^^(_;,3^ l-h-h-h,- , (56) 
T{h - hs) K2 



provided Eq. (|4q ) holds. The second contribution to (51), coming from the region p2 > pi, is simply 
evaluated by the replacements pi <-> p2 , Pi ^ P2, which give 

/ =(-1)-/ ('I'-^j), (57) 

Jp2>pi Jp2<p\ V^l ^ ^2/ 

where the parity factor stems from the change of sign of pi2 in the power with exponent —n^. 
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To derive the final expression for the Pomeron wave function in momentum space we have 
to analytically continue in the conformal weights to their physical values /13 = /i, /13 = /i and 
Ky £ TZ ■<=^ R = K*. To do this in Eq. (^) we use the relation (see Eq. (15.1.2) in [^]) 

1 



r//. . ^ 2^i(-^3, l-h;h-h3; = h{l - h)'^ 2^1(1 - h,2 - h-2;-'^) 

h3^h T{h - h-i) K2 K2 1^2 



lim 



(58) 



Putting together Eqs. ( p4| , p^ , 57,58), and rearranging some P-function factors according to the 
relation 

P(/i)P(l - /i) = (-l)"r(/i)P(l - /j) {h-h = neAf), (59) 
we obtain the final expression 



^ hil-h)Til-h)hil-h)T{l-h) ^^ 



i" (Air) 



*h-2,h-2 



-K'i ~ 2^1(1 -h,2-h;2--^) 2Fi(l -h,2-h;2--^) + (-1)"{1 ^ 2} . (60) 

K2 k\ 

As a check, we show in the following that the above function is an eigenfunction of the Casimir 
operator of the Mobius group. In the coordinate representation, in complex notation, one has 



((ri2)'9i52 + h{h - 1) 



ri2 
rir2 



0, 



(61) 



together with a similar equation in the antiholomorphic variables. In the momentum representation, 
these equations read are given by 



[{dn^ - df,^fKiK2 + Hh - 1)) E^^^{ki,K2) = 0, 



(62) 



and an analogous result holds for its antiholomorphic counterpart. From the general property of 
scaling invariance we know that 



S,;i(Ki,AC2)=K^2^I^-Xj.(l,-) 



(63) 



which is satisfied by (|6(]|). Changing the variables ki — > pi , 1^2/^1 P2 transforms (|6^ ) into 



d, 



pi 



I+P 2 
Pi 



dp, ) Pip2 + h{h - 1) ) p'l-^Ef^j^{l,p2) = 0. 



(64) 



Taking the derivative with respect to pi, we are left with a differential equation in the p2 variable 
only. We represent it in terms of a new variable y = —p2- 

(65) 



y(l - y)d'y + (2 - 2(2 - h)y)dy -{h-l)ih- 2)j E^j,il, -y) = 0. 

This is the well known hyper geometric equation. A linearly independent set of its solutions is given 
by the Kummer solutions 1*1,114 (formulas (2.9.1) and (2.9.13) in |17|): 



ui{y)= 2Fi{l-h,2-h;2-y) , u^(y) = (-y)^-^ ^F^(l - h,2 - h-2; -). (66) 

y 

These solutions match exactly the structure in (|60|), which is therefore a solution of (62). 
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One can arrive directly at ( pO[ ) by trying to construct a single-valued function from the two 
linearly independent solutions, and by further fixing the correct normalization. Note that 5-like 
terms do not appear in this approach, since we treat the holomorphic and antiholomorphic momenta 
as independent variables. The fact that, in principle, the two sectors do not simply commute (as 
one can see considering the two dimensional Poisson equation with a (^-like source) gives origin to 
the appearence of the (5-like terms. 

Let us now consider the (5- like contributions. These are present when ^ oo or /32 — > oo. 
Summing these two contributions one obtains 

El,{kuk2)=[6^'Hk,) + {-ir6^'\k,)]^2'-^^^^ q = k, + k,. (67) 

Appendix 2. Integrals appearing in the coupling to the nucleon 

In order to do the integral (|3^) we use the Feynman parametrization 
where 

D = xyA + x{l - y)B + (1 - x)C. 

In our case 

D = xy{k^ + m^) + x{l-y){{k-lf + m^) + {l-x){{k-pf + h^) = {k- x{l-y)l - {I- x)pf + R, 
where 

R = x{l- y)f + (1 - x){p^ + 6^) + xn? - (x(l - y)l + (1 - x)pf . (68) 

Shifting the integration momentum and doing the momentum integration we find the integral over 
X and y 

I = TT xdx f dy—y. (69) 



/o 'R^' 

One of the integrations (say, of y) can be done analytically. We present 

R[y) = a + Py + 'yy'^, 

where 

a = x{l - x){l - p)'^ + {1 - x)b^ + xm'^, P = x{2x - + 2x{l - x)pl, 7 = -x^/^. 
The discriminant A = 4a7 — is negative, so that the integral over y gives 

AR{1) AR{0) ''^ + 27 + ^/^/5 -^/^' ^ ' 

One finds that at x ^ the integral J behaves as so that the integration over x is 

convergent around this point. At x ^ 1 A is finite in the limit m ^ 0. However, in the limit x — > 1 
both R{0) and R{1) behave as 1/(1 — x). So we have a logarithmic divergence at x = 1, regularized 
by finite m. Evidently only the first two terms in ([70| ) lead to this divergence. Therefore we can 
safely put m = in the third (logarithmic) term. In the vicinity of x = 1 we find 
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'^^'^^ f^{l-x){p^ + h^)+m?^{l-x)[{l-pY + }fl]+m^}' ^^^^ 



The final integration over x of this term will give 



1 1 
S = I xdxJo = 77T 
l'^ 



1 A + \ 1 / (/_p)2 + 52 



In ^ 1 + Tz 77T In 



2 



(72) 



One easily checks that the singular contributions coming from /2 and /s are cancelled by the 
singular part of /i (eq. (|3^), so that the complete result is infrared finite. 

To do the numerical calculation we present J = J — Jq + Jq = Jr + Jq- The integral over x of 
the term Jr = J — Jq converges at x = 1 and can be calculated numerically. The integral over x of 
Jo is given by (|7|). 
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